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The proximity induced pair potential in a topological insulator-superconductor hybrid features 
an interesting superposition of a conventional spin-singlet component from the superconductor and 
a spin-triplet one induced by the surface state of the topological insulator. This singlet-triplet su¬ 
perposition can be altered by the presence of a magnetic field. We study the interplay between 
topological order and superconducting correlations performing a symmetry analysis of the induced 
pair potential, using Green functions techniques to theoretically describe ballistic junctions be¬ 
tween superconductors and topological insulators under magnetic fields. We relate a change in 
the conductance from a gapped profile into one with a zero-energy peak with the transition into a 
topologically nontrivial regime where the odd-frequency triplet pairing becomes the dominant com¬ 
ponent in the pair potential. The nontrivial regime, which provides a signature of odd-frequency 
triplet superconductivity, is reached for an out-of-plane effective magnetization with strength com¬ 
parable to the chemical potential of the superconductor or for an in-plane one, parallel to the 
normal-superconductor interface, with strength of the order of the superconducting gap. Strikingly, 
in the latter case, a misalignment with the interface yields an asymmetry with the energy in the 
conductance unless the total contribution of the topological surface state is considered. 

PACS numbers: 73.63.-b,74.45.+c,75.70.Tj,73.23.-b 


I. INTRODUCTION 

A three-dimensional topological insulator (3DTI) is 
an extraordinary material with an insulating band gap 
but topologically protected gapless surface states. These 
states are possible due to spin-orbit interaction and fol¬ 
low a spin-polarized Dirac-type spectrum^. Among many 
unique properties of these materials, one has triggered an 
intense interest in the condensed matter community: the 
emergence of Majorana fermions when TIs are in electri¬ 
cal contact with a superconductor^. Indeed, the topolog¬ 
ical insulator-superconductor system can be used to en¬ 
gineer an effective spinless triplet superconductor which 
is expected to host zero-energy topologically protected 
Majorana bound states^. 

Remarkably, the interplay between the spin- 
momentum locking at the surface of the 3DTI, the 
conventional pair potential from a superconductor, and 
magnetic order from either a ferromagnet or a magnetic 
field creates a very exotic induced pair potential in the 
3DTI. First, the proximity induced pair potential from 
a conventional singlet s-wave superconductor on the 
surface of a 3DTI features a superposition of singlet and 
triplet states^^^^. Owing to the symmetry of Cooper 
pairs, the singlet (triplet) components are even (odd) 
under the exchange of the spatial coordinates of the two 
electrons (even/odd in momentum). At the same time, 
they are all even under the exchange of time coordinates 
(even in frequency). Second, at the interface between 
normal and superconducting regions inversion symmetry 
is broken. Such inhomogeneity generates new pairing 
components that are odd (even) in momentum if they 
are a singlet (triplet) state of spin^’^^’^^. Therefore, 


the Pauli principle imposes that these new terms must 
be odd in frequency. Finally, as we demonstrate in 
this paper, magnetic order directly affects all these 
components and opens the possibility for the exotic 
odd-frequency terms to become dominant and produce 
strong changes in transport observables of these systems. 

Among the unconventional components of the pair po¬ 
tential, a new, very rare type of superconducting con¬ 
densate can be induced in this system: one with a spin- 
triplet component that is odd in frequency and even in 
momentum^^. Such exotic odd-frequency triplet pairing 
has been proposed to explain the long-range proximity ef¬ 
fect observed in superconductor-ferromagnet hybrids^®. 
On such systems, the odd-frequency triplet component 
can be filtered from the rest of induced components of 
the superconducting condensate since it is the only one 
resistant to disorder^®. In topological systems, spin- 
polarization can be used to detect this unconventional 
pairing. Recently, Crepin et al. have proposed a de¬ 
tection scheme based on nonlocal conductance measure¬ 
ments on the one-dimensional edge of a two-dimensional 
TI^^. Moreover, odd-frequency triplet pairing and the 
emergence of Majorana bound states are related: any 
system that supports Majorana zero-energy bound states 
also develops an odd-frequency spin-triplet component of 
the pair potentiah®^^^. 

In this paper, we consider hybrid systems between 
3DTI and superconductors like the one sketched in 
Fig. 1(a). We study the proximity effect when such sys¬ 
tems are exposed to a magnetization extended over the 
whole junction. The magnetization can be induced by 
proximity to a ferromagnetic material (b) or by an ex¬ 
ternal magnetic field (c). Previous works on 3DTI-based 


2 





Col — ^ —U) 

a a' 

k —k 

ESE 

+ 

- 

+ 

OSO 

— 

— 

— 

ETO 

+ 

+ 

— 

OTE 

— 

+ 

+ 


FIG. 1. (a) Sketch of a normal-superconductor junction on 
the surface state of a topological insulator (regions N and 
S, respectively) with an intermediate region I modeled as a 
square potential barrier (see potential profile). The angle of 
incidence 9 is measured with respect to the x-direction. (b) 
An out-of-plane magnetization is induced in the system 
by proximity from a ferromagnetic insulator, (c) An external 
magnetic field, stemming from the vector potential A (blue 
arrow) with arbitrary orientation, also generates an effective 
magnetization (red arrows), (d) Symmetry classification of 
the anomalous Green function (ai, k) . With respect to the 
exchange of the time (frequency), spin, and spatial (momen¬ 
tum) coordinates, the Green function can be even (E) or odd 
(O). The spin part can be divided into a singlet (S) or three 
triplet (T) components. The Pauli principle allows for four 
different combinations; using a “frequency/spin/momentum” 
notation: ESE, OSO, ETO, and OTE. 


normal-superconductor (NS) and Josephson (SNS) junc¬ 
tions focused only on transport signatures of such hy¬ 
brids and the conditions for the formation of Majorana 
bound states®’^°’^^“^®. Here, we study the interplay be¬ 
tween topological order and superconducting correlations 
performing a symmetry analysis of the induced pair po¬ 
tential based on the anomalous Green function. We go 
beyond previous works^’^^ by including the effect of the 
magnetization on the superconducting electrode^°. 

In the absence of magnetization, singlet and triplet 
components of the pair potential with even and 
odd frequency dependences are created at the TI- 
superconductor interface. However, the even-frequency 


singlet component is dominant and no low-energy states 
are formed^^. For finite magnetization, the odd- 
frequency triplet components can become dominant, re¬ 
sulting in a rich subgap structure with emerging low- 
energy peaks. We study their effect on the NIS conduc¬ 
tance, with I an intermediate region, and the local den¬ 
sity of states (LDOS) for two orientations of the magne¬ 
tization: in- and out-of-plane. For the latter, when the 
strength of the magnetization approaches the magnitude 
of the chemical potential in the superconductor, the odd- 
frequency triplet component becomes stronger than the 
singlet one. Consequently, a zero-energy peak emerges 
in the NIS conductance. We find the same behavior for 
in-plane case with strength comparable to the supercon¬ 
ducting gap and oriented parallel to the NIS interface. 
When the in-plane magnetization is oriented differently, 
the rotational invariance of the band structure is broken. 
As a consequence, the NIS conductance becomes asym¬ 
metric with the energy; a feature completely different 
to previous works based on ferromagnet-superconductor 
junctions on Independently of the orien¬ 

tation of the in-plane magnetization, when its strength 
is comparable to the superconducting gap, the odd- 
frequency triplet component becomes similar to or bigger 
than the singlet one. Consequently, subgap resonances 
are present in the LDOS and NIS conductance. 

This paper is organized as follows. We present the 
model for the Tl-superconductor hybrid system and de¬ 
scribe the effect of the magnetization in Sec. H. In 
Sec. HI, we describe the symmetry classification of the 
induced pair potential and study the conditions for the 
emergence of new symmetry terms for finite magnetiza¬ 
tion, which we associate with changes in the local density 
of states. In Sec. IV, we connect the previous results with 
an experimentally relevant observable like the NIS con¬ 
ductance. We present our conclusions in Sec. V. Finally, 
in the Appendix, we describe our method to compute the 
retarded Green function and the transport observables 
that we examine. 


II. MODEL 

In our setup, the surface of a 3DTI extends along the x- 
y plane. The superconducting and intermediate regions 
are located, respectively, at a: > d and 0 < x < d (see 
Fig. 1). The role of the intermediate region is to include 
any interface scattering between the normal and the su¬ 
perconducting regions. We thus assume perfectly flat 
and clean interfaces. The low-energy electron and hole 
excitations at the surface of the 3DTI are described by 
the Bogoliubov-de Gennes (BdG) equations H'^ = E'^, 
with E the excitation energy. Particle-hole symmetry 
imposes that if dli; is a solution of the BdG equations 
with excitation energy E, must also be a solu¬ 

tion. In Nambu (particle-hole) and spin space, with basis 
T' = [c-|-(k),C4,(k),c|(-k),c]^(-k)]^, with Co.(k) the an¬ 
nihilation operator for an electron of spin a =t, \, and 
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momentum k, the Hamiltonian reads as 

- ^ / h{k) - ^J.{x)ao ia2A{x) \ , . 

V -ta2Aix) ^,ix)ao-h*{-k)J ’ 

where ^J-{x) is the chemical potential and the Pauli matri¬ 
ces (Top, 2,3 act on spin space. We model the intermediate 
region as a square potential with thickness d and height 
/r/, as it is shown in Fig. 1(a). Therefore, we have 


k'N 

,a: < 0 

k-I 

, 0 < a; < d 

tJ-s 

,x > d 


We note that, due to Klein tunneling, the barrier does not 
affect the normally incident modes (with ky ~ 0)®’®. The 
pair potential in the superconducting region is A (a;) = 
A0(a; — d), with A > 0 and 0(a;) the Heaviside func¬ 
tion. Electron-like quasiparticles are described by a Dirac 
Hamiltonian as 

h(k) = vp (kxdi + ky(T 2 ) - eA- d , (2) 

with vf the Fermi velocity, e the electron charge, and A 
the vector potential (we take h = 1 and c = 1 for simplic¬ 
ity). For a constant magnetic field, the vector potential 
depends linearly on the spatial coordinates. In what fol¬ 
lows, we describe two situations where we can neglect the 
spatial dependence of A and regard it as constant. The 
effect of a uniform constant A has the same form as a 
Zeeman interaction where the vector potential plays the 
role of an exchange field. Consequently, the Hamiltonian 
in Eq. (1) is equivalent to the one used in, e.g., Ref. 23, 
when we define an effective magnetization 

= - e [{A^di + A^ds) ® ra -I- Ayd 2 0 fo] , 

where the Pauli matrices to,i, 2,3 act in Nambu space. 
The main difference with previous works^®’^® is that mag¬ 
netic order is induced in the whole junction and is not 
just limited to an intermediate region. We consider this 
a sensible approach to describe 3DTI hybrid junctions 
where proximity-induced superconductivity on the two- 
dimensional surface state can be sensitive to magnetic 
effects. 

With our choice of the reference frame (see Fig. 1), the 
momentum component parallel to the interface ky is con¬ 
served. Therefore, the solutions of the BdG equations are 
expressed as 'I'(x, 2 /) = e*^"^'I'(x,0), with sin0 = fcj^/|k|. 
A full list of the scattering states for this junction is pre¬ 
sented in Appendix A. Next, we provide the main results 
for two different orientations of the effective magnetiza¬ 
tion M. 


A. Out-of-plane effective magnetization 

We start by considering an out-of-plane effective mag¬ 
netization M = nizd^ (g) f 3 (see Fig. 1), stemming from a 


magnetic field such that A = [0,0, ^^(a;)]. Ignoring the 
spatial dependence of the vector potential is equivalent 
to assume that the magnetization is a Zeeman term with 
exchange field ruz- 

In the superconducting region, the positive branch of 
the low-energy spectrum is then given by 

£’±(k) = \J+ A"^ d- +ml± 2rz , (4) 

with rz = y/n’^Vp\'k\^ + ml {gA + A^). Although the en¬ 
ergy spectrum looks rather complicated, its analysis is 
very simple if we focus on the F point k = 0, where we 
find 

^±(0) = |m,±x/^2 + A2| . 

It is straightforward to see that both mz and the super¬ 
conducting gap A independently open a gap in an other¬ 
wise gapless Dirac spectrum. When both and A are 
finite at the same time, a competition between the su¬ 
perconducting and the magnetic gap occurs. As a result, 
a zero-energy state can be created when m\ = fjf + A^. 
This allows us to define the effective superconducting gap 

A,^ = Ay'l - (m^/l/rl)^ . (5) 

For 0 < niz < f(, the gap is real and superconductivity 
is the dominant process in the superconducting region. 
On the other hand, superconductivity is suppressed by a 
magnetic gap when > fi > A. The critical point be¬ 
tween the two regimes is thus = ^ > A. The Fermi 
energy for 3DTIs based on Hi compounds or strained 
HgTe can be estimated as /r > 50 meV. A comparable 
exchange field can be obtained by proximity from a fer¬ 
romagnetic insulator, as sketched in Fig. 1(b). In a recent 
experiment, EuS was deposited on top of the topological 
insulator Bi 2 Se 3 . Ferromagnetic order induced onto the 
surface of the TI was estimated to be around (60 — 400) 
meV/area [nm^] per applied Tesla®^. For small enough 
devices, this correction can be much larger than the in¬ 
duced superconducting gap on TIs for fields below the 
critical field of the superconductor. 

B. In-plane effective magnetization 

We now examine Eq. (3) in the heavily-doped, weak- 
field approximation with |E|,A,|m|. Since the out- 
of-plane component becomes relevant only when ^ /i, 
we can ignore it in this approximation. The resulting in¬ 
plane effective magnetization is M = mxdif 3 +myd 2 To = 
W|| ((Ti73 cos a -I- (7270 sin a), with the angle a measured 
from the kx direction [see more details in Fig. 1(c)]. The 
in-plane effective magnetization is defined from the vec¬ 
tor potential A = [Ax{z), Ay{z), 0] which we approximate 
by its value at the junction’s plane. This approximation 
is also valid for the superconducting region if we let the 
penetration length of the superconductor be much larger 
than the superconducting coherence length®^’®^. 
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In the superconducting region, the positive branch of 
the low-energy spectrum is given by 

£'±(k) Ri n • va.+ \l(fi ± nF|k|)^ , (6) 

where we have defined n = k/|k| = (cos0, sin0,0). The 
effect of the weak field is a shift in the energy, analogous 
to the Doppler shift discussed for a SQUID-like geometry 
of helical edge states in Ref. 30. Using that 

n • m = my (cos 9 cos a -|- sin 0 sin a) = TO|| cos {9 — a) , 

we can define the effective gaps 

A± = |A ± my cos (0 — a) I . (7) 

The critical value for which the effective gap closes is 
now TOy ^ A, which can be reached with rather weak 
fields. The magnetic field induces a finite Cooper pair 
(condensate) momentum which affects quasiparticles dif¬ 
ferently depending on their direction of motion. Due to 
the spin-momentum locking at the surface state of the TI, 
quasiparticles moving in opposite directions feel reversed 
magnetic fields. In other words, they travel upstream 
or downstream with respect to the superconducting con¬ 
densate. This asymmetric behavior with the direction of 
motion is not found in non-topological systems^^’^^. 


C. Retarded Green function and experimental 
observables 

We now construct the retarded Green function associ¬ 
ated to the Hamiltonian in Eq. (I). The retarded Green 
function is obtained combining all scattering states as^^ 

G^{x,x',E + iO+;9) = 

ai4'3(a;)l'f’ {x') + a 2 ^s{x)^^ {x') 

-|-Q;34'4(a:)4'f (a;')-|-a44'4(x)^|’ (x') 

/3i 4 '1 (x) (x')+/32 ^*2 (a;) (x') 

-f/33 T 1 (x) (x')+/34 2 (x) (x') 

An incoming electron (hole) from the normal region 
onto the superconductor can be ( a) Andreev reflected as a 
hole (electron); (6) normal reflected as an electron (hole); 
(c) transmitted to the superconductor as an electronlike 
(holelike) quasiparticle; and (d) as a holelike (electron¬ 
like) quasiparticle. We label this scattering process with 
the number 1 (2). The corresponding wave function is 
thus 414(2) (x). Analogously, the equivalent process where 
the incoming particle is an electronlike (holelike) excita¬ 
tion incident from the superconductor is labeled with the 
number 3 (4). In Appendix B, we show more details of 
the calculation of the Green function, the reflection am¬ 
plitudes Uj, bj, and the transmission amplitudes Cj and 
dj, with j = 1,...,4. The wave functions ^i(x) cor¬ 
respond to the conjugate scattering processes to 4'i(x). 


The coefficients and fii are determined by the conti¬ 
nuity equation obtained after integrating the BdG equa¬ 
tions around x = x', namely, 

— CTi 0 To = (9) 

VF 

G’^ix, x + 0+,E + i0+-9)- G"'(x, x-0+,E + i0+-,9) . 

The spectral density of states is then calculated from 
the retarded Green function as 

A(x, E;9) = -Im {TrG'gg(x, x,E + lO’’"; 0)} , 

where the trace is taken on the electron-electron com¬ 
ponent of the Green function in Nambu space (i.e., the 
single-particle Green function). Furthermore, the local 
density of states (LDOS) is given by 

■7r/2 

p{x,E)= J d9cos9A{x,E;9) . (10) 

-■k/2 

In the next sections, we normalize the LDOS in the su¬ 
perconductor by the LDOS calculated deep inside the 
normal region, pn{.E) = \pN + E\/vf (see Appendix G 
for more details). 

We define the conductance of the NIS junction as^®’^® 

ir/2 

Gns{E)= j A9 cos 9 [I-me, 9)\^+ \ai{E,9)\^] , 

-■n-/2 

( 11 ) 

where ai{E, 9) and bi{E, 9) are the Andreev and normal 
reflection amplitudes, respectively, for an incident elec¬ 
tron. We normalize the conductance using the constant 
value Go = Gns{E > A). 

Notice that both the conductance and the LDOS 
are averaged over all incident angles, labeled by 0 £ 
[-F/2,F/2f\ 

III. SYMMETRY OF THE INDUCED PAIR 
POTENTIAL 

Induced superconducting correlations are manifested 
in the anomalous Green function (the electron-hole com¬ 
ponent of the Green function in Nambu space), i.e., the 
2x2 matrix in spin space 

Faa'{tl,Yi-,t2,V2) = (Tt^-(ti, Ti )4'<^/(^2, r2 )) , 

with Tt the time-ordering operator and 4>cr(ii,ri) the 
fermion operator for given spin, time, and space coordi¬ 
nates. In homogeneous systems, we can define the rela¬ 
tive coordinates r = ri — r 2 and t = ti — t 2 and the Pauli 
principle imposes that = —Fa'a{—t, —r). After 

Fourier transformations, we have 


( 8 ) 

, X < x' 

, x> x' 


JAa'(^ik) = -F)tv(-w,- k) 
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We can write explicitly the spin-singlet and triplet com¬ 
ponents of the anomalous Green function in the basis of 
Eq. (1) as 

F{uj, k) = [-Fo(a;, k)(To -I- (w, k)(Tj] ia2 ■ 

Evidently, the Pauli principle must be fulfilled indepen¬ 
dently by the singlet and the three triplet components. 
With respect to the spin degree of freedom, the singlet 
component is odd and the triplet ones even. Under a sign 
change in either the frequency or the momentum, the 
functions F^{oj, k), with i' = 0,1, 2, 3, can also be classi¬ 
fied as even or odd. As a consequence, only four types of 
symmetries are allowed, namely, (i) Even-frequency spin- 
Singlet Even-parity (ESE); (ii) Even-frequency spin- 
Triplet Odd-parity (ETO); {in) Odd-frequency spin- 
Singlet Odd-parity (OSO); and {iv) Odd-frequency spin- 
Triplet Even-parity (OTE) [see Fig. 1(b) and, for more 
details, Ref. 12]. 

In our model, we assume perfectly flat and clean in¬ 
terfaces which conserve the parallel component of the 
momentum {ky for the axis choice in Fig. 1). In the 
quasi-one dimensional model, transport along these junc¬ 
tions is described by one-dimensional Green functions 
G{x,x' ,uj]ky) separately for each channel ky. For pos¬ 
itive frequencies, we construct these Green functions us¬ 
ing Eq. (8), which considers scattering solutions of the 
BdG equations with excitation energies E >Q. For neg¬ 
ative frequencies, we need to consider the time-reversal 
of the scattering processes, i.e., the E < 0 solutions of 
the BdG equations, which are provided by the advanced 
Green function. Therefore, to study the frequency de¬ 
pendence of the anomalous Green function, we combine 
the retarded and the advanced Green functions as 

F4E + isgn{E)0+]= (12) 

&{-E)E;^{E-iO+) + &{E)F^{E + iO+) , 

where if is a real variable, that coincides with the excita¬ 
tion energy for positive values, which we associate with 
the frequency, ky is parametrized using the angle of in¬ 
cidence 6 and the advanced Green function is obtained 
using G“(x, x', E — ky) = [G^{x', x,E + *0+; ky)Y . 


A. Zero-fleld analysis 

We now analyze each component of the anomalous 
Green function. We are mainly interested in local observ¬ 
ables, like the LDOS which is determined by the spatial 
variation of the symmetry of the induced pair potential. 
Therefore, we only consider the case with x = x' here. 
In Appendix B, we show that the Green function is com¬ 
posed of an edge and a bulk part, i.e., one term given by 
the scattering at the interface and another by the solu¬ 
tions far away from the interface, respectively. For each 
component of the anomalous Green function we thus have 

F, (A, 0) = (A, 0) + F,^a (A, 0) + (A, 0) . 


Fv,ao corresponds to the Green function at x = x' ^ oo, 
i.e., the bulk term inside the superconducting region. On 
the other hand, Eh and F^^b are the edge Green func¬ 
tions proportional to Andreev and normal reflection pro¬ 
cesses, respectively. In what follows, we ignore F^^b since 
it only includes rapid oscillations. 

In the absence of magnetization, for the bulk terms we 
immediately see that 


Fo,oo{E,0) 


F2,oo{E,0) 


Fi,oo{E,0) 


isgn(if) A 

2uf|cos0| y / E '^ - A2 
i sgn {E) A 

2uf|cos0| y / E '^ - A2 
El^{E,e)=0 , 


where 


Ve^ - A2 = 

f sgn (A) y/[E+i sgn {E) 0+]^ — A^ , |if| > A 
\ i sgn {E) a/A 2 — [E-l-j sgn {E) 0+]^ , [if | < A ’ 

which is an odd-function of the energy. Therefore, we find 
that Fq^oo and are even-frequency functions; they 

are thus classified as ESE and ETO terms, respectively. 
Even in the absence of magnetic order, the combination 
of a spin-orbit induced Dirac spectrum with conventional 
spin-singlet superconductivity results in a superposition 
of singlet and triplet superconducting correlations^^^^. 
The ETO term is proportional to sin0 = ky/\k\ because 
it is odd in momentum. Therefore, after averaging over 
incident angles, any odd-momentum terms must vanish. 
Consequently, the effect of ETO terms on a local ob¬ 
servable in our system is canceled by the average over 
incident angles. 

At the IS interface with x = d, the breaking of spa¬ 
tial inversion generates odd-frequency terms both for 
the singlet and triplet components^^’^^. Angle-averaging 
cancels any odd-frequency singlet term, but the odd- 
frequency triplet components remain. We thus find 


7 - 2V-E2-A2 

Fy,A{E, 0) = —a3(E;, e)G^F^^ , 

VF 

with a3{E,9) the Andreev reflection amplitude for an 
electron-like quasiparticle incident from the superconduc¬ 
tor, which results in an OTE contribution that survives 
the integration over incident angles. 

To compare the relative weight of the triplet compo¬ 
nents against the singlet one, we define the magnitude of 
the angle-averaged triplet vector 



Due to the integration over incident angles, Ft{E) repre¬ 
sents only odd-frequency triplet terms. Therefore, at zero 
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FIG. 2. (a) Evolution of p{E = 0), angle-averaged |fo|, and 
Ft with the distance inside the superconducting region from 
the interface at a: = d at zero magnetization and E = 0. (b) 
Same quantities as before as a function of for E — 0 
and X = d. The green solid and dashed lines correspond to 
the two triplet polarizations \Fi^ 2 \ oc (tt T 11)) respectively. 
(c,d) LDOS as a function of the energy for several values of 
rtiz calculated at the IS interface x = d (c), and inside the 
superconducting region, where x ^ oo (d). For all plots, we 
set pN = 10/r, pi — 20p, ps = P = 10^A, and d — Q.1 {vf/ p)- 


field, Ft{E) = \Fi,a{E)\. It is interesting to note that, in 
our spin basis, this term is given by Fi oc (ft — 11), while 
the opposite polarization, F 2 oc (ft + 11 ), vanishes. This 
is a consequence of our calculation being limited to one 
surface of the 3DTI; for the other surface the polarization 
of the Cooper pairs at the interface is reversed. 

We show in Fig. 2(a) the evolution of the angle- 
averaged quantities \Fo{E = 0)| and Ft{E = 0) as a 
function of the distance inside the superconducting re¬ 
gion (dashed blue and dotted black lines, respectively). 
The bulk ESE contribution is dominant and responsi¬ 
ble for the suppression of the density of states (red solid 
line) for x > ^. At the interface, however, we find that 
Ft{E = 0) 0. The OTE contribution vanishes inside 

the superconducting region but is responsible for the en¬ 
hancement of the LDOS close to the interface. 


hand, F^^oo is an OTE term proportional to Eruz and 
it vanishes in the absence of magnetization or for E = 0. 
Therefore, the main effect on the bulk region is deter¬ 
mined by Fq. 

The situation is different at the interface. We plot in 
Fig. 2(b) angle-averaged \Fo{E = 0)| and Ft{E = 0) 
(blue dashed and black dotted lines, respectively), as 
a function of m^, for x = d. We choose E = 0 to 
analyze the symmetry of any zero-energy peak in the 
LDOS (red solid line). The edge contributions of the 
triplet components have odd-frequency parts that remain 
finite after integrating over the incident angles. Ft(0) be¬ 
comes comparable to |Fb(0)| for a magnetization strength 
mz ~ 0.7p. Additionally, we plot the triplet components 
|^'i, 2 (A = 0)1 oc (tt T ii) using solid and dashed green 
lines, respectively. At = 0, only |F'i(0)| is present. 
For a finite niz, 1^2(0)| becomes a finite OTE term and 
increases with the field strength like jEi | at the same rate 
that the ESE term is suppressed. The LDOS at the in¬ 
terface is enhanced accordingly. With the increase of niz, 
the feature of a Majorana resonant state becomes more 
pronounced. Also, here Majorana resonances are always 
accompanied by odd-frequency pairing. 

The magnetization reduces the bulk gap A^, given in 
Eq. (5), for ruz < p. After the gap closes, the LDOS and 
the anomalous terms of the Green function are strongly 
suppressed. We show the energy dependence of the 
LDOS at the interface with x = d in Fig. 2(c) and in 
the bulk of the superconducting region with a; —>■ cx) in 
Fig. 2(d). At the interface, the LDOS features a strong 
resonance at \E\ = Az and is hnite for \E\ < Az- For 
rUz > p, the LDOS is still finite and maximum at A = 0 
[see Fig. 2(b)]. However, the value of this finite LDOS at 
the interface decreases rapidly when increasing over p 
because it is suppressed by a magnetic gap. Deep inside 
the superconducting region, where only the ESE compo¬ 
nent of the anomalous Green function is finite, we find a 
BCS-type density of states around the effective gap Az- 
For rUz > p, & magnetic gap opens destroying supercon¬ 
ductivity. 


C. Symmetry analysis for in-plane effective 
magnetization 


B. Symmetry analysis for out-of-plane effective 
magnetization 

We now study the local anomalous Green function for 
an out-of-plane effective magnetization. We present a 
list of the components of the anomalous Green function 
in Appendix D, for the perpendicular case treated here 
and for the in-plane one analyzed in the next section. 
The presence of the perpendicular magnetization affects 
all components and, in particular, induces two new bulk 
terms. Fj^oo is proportional to ky = |A:|sin0 (ETO) and 
becomes zero after integration over 9. On the other 


In the previous section, we showed that the effect of a 
magnetic held in the heavily-doped, weak-held approxi¬ 
mation is a Doppler-type shift of the excitation energy. 
This shift has an important effect on the energy bands: 
in Fig. 3(a), we plot the conduction (red) and valence 
(blue) bands in momentum space for several values of 
TO||. For reference, the zero-energy plane is marked by 
a gray area. At zero held, the spectrum has rotational 
symmetry and the superconducting gap is the same for 
jkj = kp. For hnite helds, TO|| ^ 0, the rotational sym¬ 
metry of the spectrum is broken. However, particle-hole 
symmetry is recovered using e(k) = —e(—k). The ori¬ 
entation of the NIS interfaces with respect to this bulk 
























7 



FIG. 3. (a) Energy-momentum plot of the low-energy bands 
inside the superconducting region for several values of my. (b) 
Bulk LDOS of the superconducting region of the NIS junc¬ 
tion as a function of the excitation energy E for different 
m||. (c,d) LDOS at the interface x = d with fj,i/fj, = 10 and 
d = 0.1(uf/m) S'® a function of the excitation energy E with 
different in-plane magnetizations. We show ajii = 0,0.5 for 
(c,d), respectively. For (b,c,d) we set /r = lO^A. 


FIG. 4. Plots of p{E = 0), \Fo{E = 0)|, and Ft{E = 0) 
as a function of m||. (a) Bulk results for x oo. (b,c,d) 
Results at the interface with x = d for different orientations 
of the effective magnetization: a/ir — 0,0.25,0.5 for (b,c,d), 
respectively. The solid and dashed green lines correspond to 
the polarizations tt — ij, and tt + 4-1 of th® GTE term, 
respectively. For all plots, fiN = ps = p = 10® A, pi/p = 10, 
and d = 0.1{vf/ p). 


spectrum becomes very important since the energy spec¬ 
trum can be asymmetric around ky = 0, with ky the 
momentum component parallel to the interface. Only 
the case where the in-plane field is oriented parallel to 
the normal-superconductor interface, i.e., for a = 7r/2, 
the energy bands are fc^-symmetric. Additionally, there 
is a critical magnetization my = A where both bands 
touch the zero-energy plane and the superconducting gap 
closes at some values of |k| = kp, for any orientation of 
the in-plane field. The orientation of the in-plane field 
a changes the position in momentum space where this 
closing occurs®®. 

To study the effect of the gap closing on the bulk of 
the superconducting region, we show the bulk density of 
states in Fig. 3(b) for several values of my. After av¬ 
eraging over the angle of incidence, the bulk density of 
states is independent of the orientation of the magne¬ 
tization. When 0 < my < A, the effective gaps A± 
dehned in Eq. (7) emerge and the LDOS exhibits a peak 
at E ~ I A_|-1 and is gapped for E < | A_ |. When the gap 
closes for my = A, the bulk density of states becomes 
V-shaped with two resonances at E = ±2A. For this 
case, the LDOS at the interface features a sharp zero- 
energy peak for any orientation of the in-plane field, as 
it is shown in the dashed-dotted blue lines of Fig. 3(c,d). 
Beyond this point, the zero-energy peak becomes wider 
if a < 7r/4 or splits into two otherwise. For the bulk 
LDOS, a region emerges around E = Q with a flat den¬ 
sity of states of the same magnitude as that of the normal 
region. Outside this energy range, the V-shaped profile 
becomes wider with increasing the strength of the field. 

We now study the components of the anomalous Green 
function plotting the bulk contribution of \Eo{E = 0)| 


and Ft(0) as a function of my in Fig. 4(a). There are two 
distinct regimes. For my < A, the ESE term is greatly 
enhanced, while the OTE is finite but small compared 
to it. The LDOS is suppressed which corresponds to a 
conventional gapped profile. For the critical magneti¬ 
zation my = A, all terms increase. At this point, the 
LDOS jumps from 0 to 1. From then on, the triplet vec¬ 
tor magnitude is finite and similar to the ESE term for 
my > 2A while the LDOS remains constant. Therefore, 
a finite bulk LDOS is obtained because of an equal su¬ 
perposition of OTE and ESE terms in the bulk of the 
superconducting region. 

We show in Fig. 4(b,c,d) the ESE and OTE compo¬ 
nents at the interface (a; = d), as a function of my for dif¬ 
ferent orientations of the in-plane magnetization [a = 0, 
7r/4, and 7r/2 for Figs. 4(b), 4(c), and 4(d), respectively]. 
For the symmetric orientation a = 7r/2, the OTE term 
becomes dominant over the singlet one for values equal 
or greater than the critical magnetization [Fig. 4(d)]. In 
the opposite case with a = 0 [Fig. 4(b)], the singlet term 
is always dominant over the triplet one. Independently of 
the orientation, the LDOS is greatly enhanced for values 
greater than the critical one. 

Finally, the behavior of 1 ^ 1 , 2 ] oc (ft T 4-i) strongly de¬ 
pends on the orientation of the magnetization (see green 
solid and dashed lines for 1^1,21, respectively). Indeed, 
for a = 0 we find that jFij = jF 2 j while they are very 
different for a ^ 0. Therefore, the spin polarization of 
the induced Cooper pairs can be controlled by the orien¬ 
tation of the external field®®. 
















































FIG. 5. Conductance as a function of the excitation energy 
for several effective magnetizations. For each case we include 
a sketch of the NIS junction displaying the orientation of the 
effective magnetization as a blue arrow, (a) A perpendic¬ 
ular magnetization is finite in the whole NIS junction with 
fiN ~ 10/r, fii — 20/r, and ns ~ n- The thin gray lines are 
taken at intervals O.l/r from niz = 0. (b,c) In-plane mag¬ 
netization finite in the whole junction oriented perpendicu¬ 
lar (b) or parallel (c) to the NIS interfaces. For both cases, 
Miv = MS = M m^/m — 10- (d) With the same parameters 
as (b), superposition of the conductance for both surfaces; 
sketched on the left. The black solid line corresponds to the 
my/A = 1.6 case of (b) while the dashed one is the conduc¬ 
tance of the other surface with the same parameters. In all 
plots, d = 0.1(uf/m) s-nd m = 10^A. 


IV. CONDUCTANCE SPECTROSCOPY 

Our setup, sketched in Fig. 1, is ideal for conductance 
spectroscopy measurements like the ones recently per¬ 
formed in Refs. 40 and 41. The conductance of the 
NIS junction, obtained from Eq. (11), is normalized to 
Go = Gns{.E ^ A), as it is commonly done in experi¬ 
ments. 

We start considering an out-of-plane effective mag¬ 


netization finite in the whole junction. Consequently, 
when we approach the critical point where ~ Mi 
a magnetic gap starts to develop in both normal and 
superconducting regions. In order to have propagating 
states in the normal region, we choose n-N = lOp. and 
Ms = M = lO^A. We still include a square potential 
barrier setting m/ = 20m. The results for this setup are 
shown in Fig. 5(a) for several values of rriz in the range 
0 < ruz < M- Tor ruz = 0, the tunnel conductance fea¬ 
tures a gapped profile. Due to the Klein tunneling effect, 
the conductance is not fully gapped for E < A, even in 
the presence of a strong barrier. For 0 < < Mi the 

effective gap is reduced and, as rUz approaches |J,^ the 
subgap conductance increases clearly featuring a zero- 
bias peak. The formation of this zero-bias peak coin¬ 
cides with the regime where the odd-frequency triplet Ft 
becomes dominant with respect to the ESE term |Fo| 
[see Fig. 2(b)]. For rriz > n, the gap in the excitation 
spectrum is given hy Ez = \mz — ^m^ + ^^1 and the 
conductance is zero for \E\ < Ez- 

We now consider an in-plane magnetization in the 
heavily-doped, weak-field approximation. Under this ap¬ 
proximation, if we take mtv Ms ^ l-T'l, A,my, the ef¬ 
fect of the field in the normal and intermediate regions 
is almost negligible. Assuming a setup with a symmet¬ 
ric square potential barrier with mat = Ms = M and 
M/ = 10 m, ■'^6 plot in Fig. 5(b,c) the conductance for sev¬ 
eral values of my. We show an in-plane magnetization 
perpendicular to the NIS interfaces (a = 0) in Fig. 5(b) 
and parallel to them (a = 7r/2) in Fig. 5(c). In both 
cases, the effect of the field is to split the superconduct¬ 
ing gap into two effective gaps A± defined in Eq. (7). 

For a = 7r/2, the resonance at U = A for my = 0 
remains pinned to the smallest of the gaps, but the one 
at the biggest gap is smeared by the angle average. Due 
to Klein tunneling, these resonances can not be clearly 
resolved in the conductance, although they appear in the 
LDOS, as shown in Fig. 3(c,d). The smallest of the gaps 
closes for my > A and a zero-energy peak appears in the 
conductance. It is in this regime where the odd-frequency 
triplet term is greater than the singlet one, as shown in 
Fig. 4(d). 

Strikingly, for a ^ 12, the conductance is asymmetric 
with the energy. This asymmetry is a consequence of the 
distortion of the energy bands by the magnetic field, as 
shown in Fig. 3(a). The field breaks the rotation sym¬ 
metry of the energy bands and the modified bands have 
a point-like symmetry that is determined by the orienta¬ 
tion of the in-plane magnetization. Since we study trans¬ 
port in the direction perpendicular to the NIS interfaces, 
when the magnetization is aligned with the interface, i.e., 
a = 7r/2, the symmetry is recovered [see Fig. 5(c)], while 
it is lost for any other orientation. 

Owing to the special particle-hole symmetry of our sys¬ 
tem, we find that, for a ^ 7r/2, Andreev reflection is 
very asymmetric with the energy while normal reflection 
is symmetric. This effect holds independently of the bar¬ 
rier strength provided that the effective magnetization is 
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misaligned with the NIS interface. Thus, as an exam¬ 
ple, we can consider a junction with fxj = jis fJ-N- 
For such junctions, normal and Andreev reflections take 
place separately at each side of the intermediate region. 
In the setup of Fig. 1(a), when \E\ < A, normal (An¬ 
dreev) reflection is the only scattering process at x = 0 
{x = d). The reflection amplitudes are then 


ai{E,e) 

bi{E,9) 

t{E,e) 


2 cos cos 9 s 

t(Ej) 

.sin^AT — sin 0s 


r- , 
(i-r+r_) 


t[E,9) 

1 -I- cos {9n + 9s) 

-r+r_ [1-cos(0w-0s)] 


(14a) 

(14b) 


with 0Ar_s = sin ^{ky/fiN,s) and 


r± = 


e± + - |AP 

e± = E ± mil cos(0s — <a) 


Due to the spin-momentum locking at the surface of 
the 3DTI, the magnetic field affects differently quasipar¬ 
ticles moving in opposite directions^®. That is the reason 
why Andreev reflection is asymmetric with the magne¬ 
tization, being proportional to F_ while normal reflec¬ 
tion is not, since it goes with (1 — F+F-). As a con¬ 
sequence, the conductance becomes more asymmetric as 
we increase my, as shown in Fig. 5(b). The energy range 
where Andreev reflection is enhanced also becomes wider. 
In the hgure, this happens always for positive energies. 

Finally, the asymmetry in the conductance is also 
partly a consequence of having limited our analysis of 
transport to only one surface of the 3DTI. In Eq. (1), the 
Dirac Hamiltonian chosen to describe electron-like quasi¬ 
particles is /i(k) = [up’k—eA]-(T. With this choice, quasi¬ 
particles have positive helicity, i.e., their spin points in 
the direction of their momentum. For the opposite sur¬ 
face state of the 3DTI, however, electron-like quasiparti¬ 
cles are described by /i(k) = [— n_Fk — eA] • a. Therefore, 
the sign of their helicity changes and their spin points 
in the opposite direction of their momentum. As a con¬ 
sequence, the effect of the in-plane magnetization is re¬ 
versed, with r± —>■ Fzp. The normal reflection amplitude 
is unchanged, but the Andreev reflection amplitude be¬ 
comes proportional to F+ in Eq. (14). The asymmetry of 
the conductance is reversed, as it is shown by the black 
solid and dashed lines of Fig. 5(d), corresponding to the 
TO||/A = 1.6 case of Fig. 5(b). If we consider a conduc¬ 
tance measurement that includes opposite surfaces of the 
3DTI, the symmetry with the energy is recovered. For 
the same parameters as Fig. 5(b), we show such conduc¬ 
tance spectroscopy in Fig. 5(d). 


V. CONCLUSIONS 

We have analyzed the superconducting proximity ef¬ 
fect at 3DTI-superconductor hybrids in the presence of 
an effective magnetization. We find that the interplay be¬ 
tween the magnetization and the spin-momentum locking 
of the surface state gives rise to interesting odd-frequency 
triplet terms in the induced pair potential. Contrary to 
the usual even-frequency triplet terms, averaging over the 
angle of incidence does not cancel these components, and 
their effect can be observed in experimental transport ob¬ 
servables such as the NIS conductance and LDOS. 

For an out-of-plane magnetization, the OTE terms 
increase with the strength of the field in the regime 
/i > niz > 0. At the same time, the singlet term of the 
induced pair potential is suppressed. As a consequence, 
the conductance evolves from a conventional gapped pro¬ 
file for niz = 0 into a zero-bias conductance peak when 
Such strong magnetization can be induced in 
the system when a ferromagnetic insulator such as EuS 
is deposited on top of the NIS junction®^, in a setup like 
the one sketched in Fig. 1(b). 

If we consider the effect of an orbital magnetic field on 
the junction, we find a similar gap to zero-energy peak 
transition on the conductance for in-plane magnetization 
parallel to the NIS interface. The main difference is that 
the critical value that determines the emergence of a zero- 
energy peak is now rn|| = A, where the energy bands 
feature an indirect closing of the gap. For my > A, the 
OTE term is greater than the ESE one at the interface 
and similar in magnitude in the bulk superconducting 
region. As a result, we have connected the emergence 
of odd-frequency triplet in the induced pair potential on 
the surface of the 3DTI with a distinctive profile of the 
conductance. 

When the in-plane magnetization is not oriented par¬ 
allel to the interface, the conductance becomes asym¬ 
metric with the energy. This is a manifestation of the 
breaking of the rotational symmetry of the band struc¬ 
ture of the combined 3DTI-superconductor system. Only 
a transport measurement that probes one surface state of 
the 3DTI can detect the asymmetry of the conductance. 
Such setup, would be a perfect platform to probe the 
coexistence of topological order and superconductivity. 

Our results show that 3DTI-based NIS junctions fea¬ 
ture a very rich induced superconducting pair potential 
that displays the elusive odd-frequency triplet compo¬ 
nent. Using very basic ingredients such as a conventional 
superconductor and an external magnetic field or fer¬ 
romagnetic insulator, standard experimental techniques 
like conductance spectroscopy can be used to detect sig¬ 
natures of OTE superconductivity. 
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Appendix A: Scattering states 

In this appendix, we present the solutions of the BdG 
equations 774' = £'4', for the two orientations of the 
effective magnetization considered in the main text. In 
the normal regions with A = 0, the general solution of 
the BdG equations is 


^(x <d)= A+V'+feeS*"'" 


Analogously, for the superconducting region, the wave 
vectors are given by kj = — {vpkyY/{vp): 

with j = 1,2, 


n(2)± 


J- (—± ErUz 
fi±mz 


_ f sgn (£ ) v ^£ 2 _ A2 ,|£| > A , 

VA2 - £2 ,|£|<A, 


and A^ = A-^/T — (m-zZ/i)^. The corresponding wave 
functions in Eq. (A2) are 


V'ifc, = [ri,±£i±Fi,TC£^i±,I]^ , (A4a) 

= [r2,±£2±F2,TCi^2±,I]^ , (A4b) 


with 


£ 


1 ( 2 )± 


ri.2 = 




VF(ki(2) T iky) 

/i£ ± fiflz 

A(fi - ruz) 

^i + ruz 


C = rir2 = 


fj.-mz 


where the coefficients [Bjj) represent the amplitudes 
for electrons (holes) propagating to the right and left, 
respectively. 

On the other hand, for the superconducting region we 
find 

^(x >d)= A+V'+fc^e^'^^" + As^_fc,e-*'='" 

, (A2) 

where the coefficients A^ (£^) now label the amplitudes 
for right and left moving electron-like (holes-like) quasi¬ 
particles. 


1. Junctions with out-of-plane effective 
magnetization 

In the normal state regions with A = 0, the wave vec¬ 
tor is given by ke,h = \/Ze+Zg- - {vpky)^/{vp), with 
Zg± = /i + £ ± TTiz and Zh± = — E E ruz . The wave 

functions in Eq. (Al) are then given by 

= [l,±e±*®G0,0]^ , (A3a) 

V'ifc. = [0,0,1, , (A3b) 

where we have defined the phase factors 

±20= i ±20;, ^ Sh'lky 

9~ ^ _ y -^ g=n2t7;i -y 

^ IZe+l ’ |£,-| ’ 


2. Junctions with in-plane effective magnetization 

In the normal region, where A = 0, the wave vector 
component perpendicular to the interface becomes = 

ke,h±rna;/vF with vpke^h = ^(/i ± Ef-{vFky ± m^)^. 
The superscript > (<) labels right (left) movers along the 
x-direction. Electrons and holes have an opposite shift 
by the field. When they are coupled by the superconduc¬ 
tor, this shift becomes very important. Therefore, the 
component of the in-plane magnetization perpendicular 
to the interface, nix = wy cos a, discriminates between 
particles that move with or against the stream of Gooper 
pairs. The general solution of the BdG equations for the 
normal state region is given in Eq. (Al) with the same 
wave functions defined in Eq. (A3), but the new phase 
factors 

^ vpke^h + i {vpky ± niy) 

\d±E\ 

In the heavily-doped regime :s> £, A, we have 9 = 


sin ^ {ky/^), ki ^2 ^ kFx ± *^ 2 , and —ki ^2 ^ 
with 

Cl AP — (£ ± n • m)^ 

■' kpx ^ 5 

vfcos6 

5 

and kFx = cos0. Under this approximation, the wave 

functions in Eq. (A2) are 

= [l,±e±*^T^^e±*^^^]^ 

, (A5a) 

= [F^,±F^e±*^Te±*^l]^ 

(A5b) 


and Sh = sgn(/r — £). 
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with 

P ^_ A* 

* (ii ± n • m) + yJ{E ±n- m)^ ~ I Ap 


Appendix B: Green function techniques for Dirac 
systems 


For an incoming electron from the normal region, the 
wave function is 


'Si{x)=< 


a: < 0 

(Bl) 



X > d 


For the other processes, we find 



^'2(a;)=i 

' e *'='*^V'-fcfe+a 2 e 

, X < 0 

1 

(B2) 

1 

C2e“*'=2^'!/'-fe2+t^2e*'^'i^V^+fci , 

X > d 


1 


X < 0 


^'3(x)=< 


, X > d ’ 

(B3) 

1 



and. 


X < 0 


4'4(x)='! 

i +64e-*'=^"V'-fe2 

X > d 

(B4) 

1 


At the intermediate region, we use the wave functions 
d/,(0 <x<d)= (B5) 

which correspond to the normal region solutions with the 
change p ^ pi and where j = 1,... ,4 labels the pro¬ 
cesses. 

The corresponding reflection and transmission ampli¬ 
tudes are obtained inserting Eqs. (A3), (A4), and (A5) 
into the boundary conditions 

vI/,(0-) = vI/^(0+) , vI;^(d-0+) = vI/,(d+0+) . (B6) 

Next, we insert the resulting wave functions for the 
scattering processes in the retarded Green function de¬ 
fined in Eq. (8). The wave functions correspond to 

the conjugate scattering processes to and are solu¬ 

tions of H*(k) with the change k — ?> —k. The conjugated 
wave vectors are obtained by the transformation 

'^±kAky) = : 

which is a parity transformation. 

We now consider the Green function of the normal re¬ 
gion, imposing that ke^h > 0. We substitute the solutions 
with a; < 0 of Eqs. (Bl), (B2), (B3), and (B4) into Eq. (8) 
and apply the continuity condition of Eq. (9) to obtain 


X 


Ake\x-x'\ 


Ake.\x-x'\ 
COS de. 


'I’-tJ’-t. + I’l ' 




COS 0e 

-i(kex'-k^x) J. 

«2- 

-i(kcx'-k^x) J. 

— i’+kkA+k 


Or 


-i(kex-kf^x') 


-ikklx-x'\ 


COS Gf, 


'^-ke'^-kk 

Akh(x + x') 


'+ke 


COS 9h 


\ I , t , , Q'ikhix + x ) I 

-^-kkA-kk + ^2 i’+kkA-kk 

-i(kex-k^x') J. 

« 2 — 

ikk\x-x'\ , , .]■ , ^ikk{x + x') , I 


COS Qh 


I / , t I L + ^ ) , /1 

+b2 ^+kki’-kk 



(B7) 


Finally, the Green function on the superconducting side of the junction, with Re{A:i_ 2 } > 0, is 


GAx^x') = 


Ai 



-\~A.2 


Bl 



+B 2 

< 



(B8) 


where 


Ai =Si 


» ^1- tA 

VF ki Fi — F 2 

i Z 2 - TA 

Up ^2 Ti— F2 


for an out-of-plane magnetization and 


Ai = B 2 
A 2 = Bl 


i 1 

vp cos 9 \ — T\ 

i 1 

UFCOS0 1 —r?_ 


A 2 — B 2 
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for the in-plane magnetization in the heavily-doped, 
weak-field approximation ^ A,i5,m||) and with |0| < 
7 r/ 2 . For the Green function with Re{fci 2 } < 0, the 
definition of the coefficients Ai ^2 and B 12 changes sign. 
Consequently, the sign change in 1 /A:i _2 or, equivalently, 
1/ COS0 with 7 r /2 <0 < 37 r/ 2 , is canceled. 

The Green function can be trivially separated into a 
bulk contribution, defined far away from the interface, 
and an edge term which contains the scattering at the 
interface. This is, G'^{x,x') = G^{x,x') + G\{x,x') + 
Gg(a;,x'), where we have divided the edge term into 
the part that is given by Andreev reflection processes 
\G\{x,x')\ and the one that is given by normal reflec¬ 
tions [G^b{x,x')]. 


Appendix C: Local density of states 


superconducting region comparable to the superconduct¬ 
ing coherence length. 

When we consider an in-plane magnetization, the pre¬ 
vious result is changed to 


7r/2 


p{E, X > d) = - [ 

nvF J 

i sgn {E) 


A6 cos 9 


-7r/2 


Im 

r+ 

i-n 


cos 9 


1 - r^r 2 


03 


r_ 

•04 


i-r 2 


e^f_b 3 2ikix _ T— 


1 -r 


1 -r 


+ 

^3 p2ikix 

" 2 ” ^ 


e “"■^■“64 ,E>0 
i^e*®e-2*fc2^64 , A < 0 


(C4) 


tan 0 


The electronic LDOS is obtained from the retarded 
Green function using Eq. (10). In the normal region is 
thus given by 


Pn{E) = p{E,x ->• - 00 ) 

Tr/2 

= - [ d 0 cos 01 m|—-^rr i (Cl) 

TTVF J L COS 9e J 

-■n-/2 

for the bulk and 

7r/2 

Pb(E,x) = -/ d0cos0 

TTVF J 

Tr/2 

xIm{6i(A,6l)tan6lee"*®-=e"2*'^'^} , (G2) 

for the edge. The LDOS in the normal region is thus 
piE,x < 0) = pn{E) + pb{E,x). 

Analogously, in the superconducting region we find, in 
the absence of magnetization. 


t/2 


p{E, X > d) = -- 


d0 cos 9 


ttvf 


(C3) 


Im 


i sgn (E) /I + r 


COS0 


--n/2 

2 r 


i-r 2 i-r 2 


oae 


— 2k,x 


+ 


tan 0 

T^T 2 


iS 2ikixi 




2^—i6^—2ik2X 


hi) 


We have neglected the rapidly oscillating terms propor¬ 
tional to the normal reflection amplitudes h'i^i{E,9) in 
the edge contribution. The bulk contribution of Eq. (C3) 
reduces to the BGS density of states p{E,m = 0) = 
Re { — A^}. We plot in the left panel of Fig. 6 (a) 

the LDOS for zero field as a function of the energy and 
the distance from the interface inside the superconduct¬ 
ing region. The LDOS is finite for |i5| < A at the inter¬ 
face {x = d) and decays to zero for distances inside the 


We plot in Fig. 6 (b,c) the LDOS as a function of the 
energy and the distance from the IS interface at x = d 
for several values of m|| and for a = 0 and a = 7 r/ 2 , re¬ 
spectively. The effect of the in-plane magnetization is to 
split the superconducting gap into two, A —>■ A 4 -. When 
my < A, the LDOS at the interface clearly shows four 
resonances at \E\ = Aj- which become, inside the super¬ 
conducting region, fully gapped for \E\ < A_ with sharp 
resonances at \E\ = A+ [see left panels of Fig. 6 (b,c)]. 
For my = A, A_ = 0 and the LDOS inside the super¬ 
conducting region adopts a V-shaped profile. At the in¬ 
terface, a peak at A = 0 appears which has a long-range 
decay if a = 0 (|a: — d| ^ 10 ?)- For a magnetization ori¬ 
ented parallel to the interface (a = 7r/2), the peak decays 
as \x — d\ ^ This peak disappears when my > A if 
a > 7r/4, but it is still present and displays a long-range 
decay otherwise. 

Finally, for an out-of-plane magnetization, the LDOS 
is given by 

p{E,x > d) = - [ dfci, Im {sgn (A) (G5) 

TTVF J 

X [AiF?(l + E 1 +E 1 -) + A 2 F 2(1 + E 2 +E 2 -) 
+riF 2 (Aia 2 -l- A2a4)e~^''*^ 

J 1 + Ex+E2- , E > 0 1 

\ 1 + E^_E2+ ,A < 0 j 

-A2TI{1-El_)e-'^^^^^hi]] . 

We show the LDOS results for a perpendicular magneti¬ 
zation close to the closing of the gap in the central and 
right panels of Fig. 6 (a). Before the closing of the gap 
(central panel with m^ = 0.9), the LDOS shows a reso¬ 
nance at the energies corresponding to the effective gap 
A;; and is greatly enhanced at zero energy close to the 
interface. This enhancement of the LDOS, however, de¬ 
cays fast inside the superconducting region and disap¬ 
pears at a distance comparable to the superconducting 
coherence length. After closing the gap (right panel with 
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E/A 


mu/A = 


0.4 


E/A 


E/A 


m\\/A = 1.0 


E/A 


1.6 0 = 7^2 


V 


E/A 


m^//i = 0.9 


FIG. 6. LDOS inside the superconducting region normalized to the bulk density of states in the normal region pjv as a function 
of both the excitation energy and the distance from the interface {x = d). (a) Out-of-plane magnetization. From left to right, 
= 0,0.9,1.1. For the three panels, pjv = 10/r and p/ = 20p. (b) In-plane magnetization perpendicular to the NIS 
interface {a = 0). From left to right, m||/A = 0.4,1.0,1.6. (c) In-plane magnetization parallel to the NIS interface {a = 7r/2). 
From left to right, m||/A = 0.4,1.0,1.6. For (b) and (c), pjv = /r and p/ = 10/r. For all plots, p = lO^A and d — 0.1{vF/fJ-)- 


niz = 1 . 1 ), superconductivity is strongly suppressed and 
only a zero-energy peak on the LDOS survives at the 
interface. This peak decays inside the superconducting 
region within 10 “^^. 


In the absence of magnetization, the Andreev reflection 
probabilities adopt the simple form 03 = 04 = da = 04 = 
—r and the normal reflections are da = 64 = 63 = 64 = 0 . 
The components of the anomalons Green fnnction are 


Appendix D: Anomalous Green function 

We now analyze the electron-hole component in 
Nambu space of Eq. (B 8 ). For simplicity, we only con¬ 
sider NS jnnctions with no intermediate region. More¬ 
over, we only show the case with Rejfca,} > 0. The 
symmetry classification remains the same for more com¬ 
plicated NIS junctions. 
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thus given by 


Fo{E,e) = 


i sgn [E) 


2vf cos 6 
F3{E,0)=O, 
E^^{E,e)=EuiE,e) 
sgn (E) 


2vf cos 6 
id 


2 r 


i_r 2 i-r 2 


2 r 


2vf 


Use 


i_r 2 

— 2kx 


1 + r" 

036 


1 _i_ r 2 

-L ^ i- 


sin( 


with a = +1,-1 for and Fn, respectively. We 
have defined F = A/(£' + r2), k = —iil,/{vFC0s9), and 
il = — A^. The singlet term Fo{E,9) is even in 

frequency and spatial dependence; therefore, it is classi¬ 
fied as ESE. The first term of F^„{E,9), with a =t, 4'5 is 
equal to Fq{E, 9) multiplied by sin 9, which makes it even 
in energy but odd in spatial dependence, hence classified 
as ETO. The second term of F„„{E,9) is odd in energy 
and even in spatial dependence and is classified as OTE. 

In the main text, we use the triplet components Ei ,2 = 
F^^ =F F^^. Due to the sign change for the edge part 
proportional to the Andreev reflection amplitude 03 , F 2 
is equal to the bulk ETO term multiplied by sin 9 while 
Fi is given by the edge OTE part. After averaging over 
the angle of incidence, only the ESE and OTE terms are 
non-zero and we obtain the behavior shown in Fig. 2(a) 
of the main text. 

We now consider an in-plane magnetization. For NS 
junctions with no intermediate region, we find 03 = 0,4 = 
—r+, a 4 = 5-3 = —r_, and 63 = 64 = 63 = 64 = 0. The 
singlet component is 


Fo{E,9) = 


i sgn {E) 
2vf cos 9 
r" 


(t 

(t 


1-n 


le 


— (ki—K2)x 


tO-s 


1 -r 

r 


^04 


^03 


i-r 


^0-4 


1 — r?_y ' 2vf cos 9 
E>0 
E <0 


(Dl) 


which reduces to the previous result for my = 0 , where 
r_|_ = r_ = r, and is still classified as ESE. For the 
triplet components, we find that F3{E,9) = 0 and 


Faa{E,9)= 


r_ 


sgn (A) / r+ 

2cose U-r+ ' 1-r?. 

+ ^ i_r^ + i_rii , E > 0 

— ^ i-T^ *^3 F , A < 0 


^-{k,i-K2)x 

2vf cos 9 


sine 


2vf \1-Tl 1 - TlJ 

I + (i^^3 - 1 ^ 04 ) 




2vf 


, E>0 
, E<0 


(D2) 


As before, the term proportional to sin 9 is even in energy, 
thus belonging to ETO classification, while the term pro¬ 
portional to a is odd in energy and is classified as OTE. 

Finally, we consider the out-of-plane magnetization. In 
this case, the reflection amplitudes adopt a rather com¬ 
plicated form. For simplicity, in the following analysis we 
only consider the terms proportional to Andreev reflec¬ 
tion amplitudes, which fulfill ai = di. The terms coming 
from normal reflections only add some rapid spatial os¬ 
cillations. Under these approximations, the spin-singlet 
component of the anomalous Green function is 


Fo{E,ky) = 


i sgn {E) 


4?;F(ri-r2)V ki 

-2^^x /z,.r 


^1+ F C-^l- ^ -^2-1- + (Z 2 - 


le 


ivF (Fi - F2) V Tifei 
X sgn (E) 


1-1-4 2 , ^-^2-1^1 

03 + C 7- CI 4 


^2k2 


T , kik2 ki-k2 

J- ^ ^ 


■Z^l-l-^ 2 - 


"^^ 1 +^ 2 - 


(D3) 


The bulk part of the singlet component is even in 
energy and belongs to ESE classification. For the 
edge part, proportional to exp(— 2 K 2 a;), with = 

iF — ml — Ki —i(ki — k 2 )/ 2 , we find another 

ESE term together with an OSO term. The latter, pro¬ 
portional to ky, is only present when there is electron-hole 
asymmetry and vanishes if /i 3 > A, A. 

When the out-of-plane magnetization is finite, we find 
a new triplet component, which was zero in the previous 
analysis, namely. 


F3{E,ky) = - 


i sgn (E) 


Em-z 


2vf (El — F2) fi — m 

■e- 2 «.. /Z1+F2 , Z2-T 

-as + C 


1 

k^ 


1 


4i;f (Fi - F 2 ) V Eifci 
kik2 


sgn (A) 1 - 


^ 14 -^ 2 - 


— ik-,- 


F 2 fc 2 

ki — k2 


-04 


^ 1 -f ^ 2 - 


(D4) 


The bulk part, which classifies as OTE, is zero for E = 0. 
The edge part has components from both ETO and OTE. 
For the other triplet components, we find 


FniF,ky)= - ky sgn {E)C 


ki + k2 


2vf (Fi — F2) 


-b 


i(e 


4vf (El — E 2 ) 
+iky sgn (E) 


E 2 El 

^03 - 4^04 

4 1 4 2 


E 2 ^ - El 

+ CtTI—«4 


and 


FuiF,ky) = -ky sgn {E) 


Eifci E 2 fc 2 

ki + k2 


2vf (El — E2) 


le 


Avf (El — E 2 ) 
+iky sgn (A) 


Zi+ T2k2 CZ 2 - Eifci 

-®3- , 

Zi-i- i 2 fc 2 


(D5) 


(D 6 ) 


C^ 2 - Eifci 


04 


Z\F E 2 C-^2- El 

■03- 


C^ 2 - Eifci 


h+ ^2k2 


As described in the main text, both bulk terms have ETO 
symmetry and are canceled after averaging over incident 
angles. 







































































15 


^ X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 

( 2011 ). 

^ J. Alicea, Reports on Progress in Physics 75, 076501 
(2012); C. Beenakker, Annual Review of Condensed Mat¬ 
ter Physics 4, 113 (2013). 

® L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 (2008); 
J. D. Sau, R. M. Lutchyn, S. Tewari, and S. Das Sarma, 
ibid. 104, 040502 (2010); R. M. Lutchyn, J. D. Sau, and 
S. Das Sarnia, ibid. 105, 077001 (2010); J. Alicea, Phys. 
Rev. B 81, 125318 (2010); Y. Oreg, G. Refael, and F. von 
Oppen, Phys. Rev. Lett. 105, 177002 (2010); A. C. Potter 
and P. A. Lee, ibid. 105, 227003 (2010); Phys. Rev. B 83, 
184520 (2011); E. Prada, P. San-Jose, and R. Aguado, 86, 
180503 (2012); M. Houzet, J. S. Meyer, D. M. Badiane, 
and L. 1. Glazman, Phys. Rev. Lett. Ill, 046401 (2013); 

R. P. Tiwari, U. Ziilicke, and C. Bruder, 110, 186805 
(2013); New Journal of Physics 16, 025004 (2014); S. Rex 
and A. Sudbp, Phys. Rev. B 90, 115429 (2014). 

N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 

® L. P. Gor’kov and E. 1. Rashba, Phys. Rev. Lett. 87, 
037004 (2001). 

® A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Lud¬ 
wig, Phys. Rev. B 78, 195125 (2008). 

^ T. Yokoyama, Phys. Rev. B 86, 075410 (2012). 

® G. Tkachov and E. M. Hankiewicz, Phys. Status Solidi B 
250, 215 (2013). 

® G. Tkachov, Phys. Rev. B 87, 245422 (2013). 

G. Tkachov and E. M. Hankiewicz, Phys. Rev. B 88, 
075401 (2013). 

P. Burset, F. Keidel, Y. Tanaka, N. Nagaosa, and 
B. Trauzettel, Phys. Rev. B 90, 085438 (2014). 

12 Y, Tanaka and A. A. Golubov, Phys. Rev. Lett. 98, 037003 
(2007); M. Eschrig, T. Lofwander, T. Champel, J. Cuevas, 
J. Kopu, and G. Schon, Journal of Low Temperature 
Physics 147, 457 (2007); Y. Tanaka, A. A. Golubov, 

S. Kashiwaya, and M. Ueda, Phys. Rev. Lett. 99, 037005 
(2007); Y. Tanaka, Y. Tanuma, and A. A. Golubov, Phys. 
Rev. B 76, 054522 (2007); Y. Tanaka, M. Sato, and N. Na¬ 
gaosa, Journal of the Physical Society of Japan 81, 011013 
( 2012 ). 

A. M. Black-Schaffer and A. V. Balatsky, Phys. Rev. B 86, 
144506 (2012); 87, 220506 (2013); 88, 104514 (2013). 

V. L. Berezinskii, JETP Letters 20, 287 (1974); Pisma Zh. 
Eksp. Teor. Fiz. 20, 628 (1974); A. Balatsky and E. Abra¬ 
hams, Phys. Rev. B 45, 13125 (1992). 

F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Phys. Rev. 
Lett. 86, 4096 (2001); T. S. Khaire, M. A. Khasawneh, 

W. P. Pratt, and N. O. Birge, 104, 137002 (2010). 

For more information about odd-frequency triplet pairing 
on ferromagnet-superconductor hybrids see F. S. Bergeret, 
A. F. Volkov, and K. B. Efetov, Rev. Mod. Phys. 77, 1321 
(2005); M. Eschrig, Phys. Today 64, 43 (2010). 

F. Crepin, P. Burset, and B. Trauzettel, Phys. Rev. B 92, 
100507 (2015). 

Y. Asano and Y. Tanaka, Phys. Rev. B 87, 104513 (2013). 
V. Stanev and V. Galitski, Phys. Rev. B 89, 174521 (2014). 

X. Liu, J. D. Sau, and S. Das Sarma, Phys. Rev. B 92, 
014513 (2015). 

H. Ebisu, K. Yada, H. Kasai, and Y. Tanaka, Phys. Rev. 


B 91, 054518 (2015). 

M. Snelder, A. A. Golubov, Y. Asano, and A. Brinkman, 
Journal of Physics; Condensed Matter 27, 315701 (2015). 
23 Y, Tanaka, T. Yokoyama, and N. Nagaosa, Phys. 
Rev. Lett. 103, 107002 (2009); J. Linder, Y. Tanaka, 
T. Yokoyama, A. Sudbp, and N. Nagaosa, Phys. Rev. B 
81, 184525 (2010); A. Soori, O. Deb, K. Sengupta, and 
D. Sen, 87, 245435 (2013). 

J. Linder, Y. Tanaka, T. Yokoyama, A. Sudbp, and N. Na¬ 
gaosa, Phys. Rev. Lett. 104, 067001 (2010). 

C. T. Olund and E. Zhao, Phys. Rev. B 86, 214515 (2012). 
M. Snelder, M. Veldhorst, A. A. Golubov, and 

A. Brinkman, Phys. Rev. B 87, 104507 (2013). 

J. Nussbaum, T. L. Schmidt, C. Bruder, and R. P. Tiwari, 
Phys. Rev. B 90, 045413 (2014). 

B. Lu, P. Burset, K. Yada, and Y. Tanaka, Superconductor 
Science and Technology 28, 105001 (2015). 

B. Lu, K. Yada, A. A. Golubov, and Y. Tanaka, Phys. 
Rev. B 92, 100503 (2015). 

G. Tkachov, P. Burset, B. Trauzettel, and E. M. Han¬ 
kiewicz, Phys. Rev. B 92, 045408 (2015). 

P. Wei, F. Katmis, B. A. Assaf, H. Steinberg, P. Jarillo- 
Herrero, D. Heiman, and J. S. Moodera, Phys. Rev. Lett. 
110, 186807 (2013). 

M. Fogelstrom, D. Rainer, and J. A. Sauls, Phys. Rev. 
Lett. 79, 281 (1997); Y. Tanaka, Y. Tanuma, K. Kuroki, 
and S. Kashiwaya, Journal of the Physical Society of Japan 
71, 2102 (2002); Y. Tanaka, T. Yokoyama, A. V. Balatsky, 
and N. Nagaosa, Phys. Rev. B 79, 060505 (2009). 

P. Burset, W. Herrera, and A. Levy Yeyati, Phys. Rev. B 
80, 041402 (2009). 

W. L. McMillan, Phys. Rev. 175, 559 (1968); A. Fu¬ 
rusaki and M. Tsukada, Solid State Communications 78, 
299 (1991); Y. Tanaka and S. Kashiwaya, Phys. Rev. B 
53, 9371 (1996); S. Kashiwaya and Y. Tanaka, Reports 
on Progress in Physics 63, 1641 (2000); W. J. Herrera, 
P. Burset, and A. L. Yeyati, Journal of Physics: Con¬ 
densed Matter 22, 275304 (2010). 

G. E. Blonder, M. Tinkham, and T. M. Klapwijk, Phys. 
Rev. B 25, 4515 (1982). 

36 Y, Tanaka and S. Kashiwaya, Phys. Rev. Lett. 74, 3451 
(1995). 

The one-dimensional edge of a 2DTI, which is equivalent 
to the 6 — 0 case considered here, presents a well-defined 
spin polarization. For the 2D surface state, the modes with 
6^0 usually spoil this neat effect. 

The gap closing can greatly affect the edge states of 
the two-dimensional TI, as recently described in R. W. 
Reinthaler, G. Tkachov, and E. M. Hankiewicz, Phys. Rev. 
B 92, 161303 (2015). 

The orientation of the in-plane vector potential A and, 
hence, that of the effective magnetization m||, is shifted 
from that of the magnetic field by — 7r/2. 

A. D. K. Finck, C. Kurter, Y. S. Hor, and D. J. Van Har¬ 
lingen, Phys. Rev. X 4, 041022 (2014). 

M. Snelder, M. Stehno, A. Golubov, C. Molenaar, 

T. Scholten, D. Wu, Y. Huang, W. van der Wiel, 
M. Golden, and A. Brinkman, (2015), arXiv:1506.05923. 



